Myc controls NK cell development, IL-15-driven expansion, and translational machinery

Myc ablation in vivo impacts immature NK cells’ ribosomagenesis and development. Accordingly, mice with NK cells lacking Myc exhibit impaired anticancer immunity.


Introduction
NK cells, which are important to control selected cancers and infections, express multiple germ-line-encoded activating and inhibitory receptors. The latter are engaged by major histocompatibility complex class I (MHCI) molecules, whose expression can be lost by hazardous cells trying to evade T cell responses. Conversely, activating receptors are engaged by ligands expressed by stressed cells, including infected and transformed ones. NK cell cytotoxic activity therefore depends on the imbalance towards activating signals (Narni-Mancinelli et al, 2011). Recently, the generation of chimeric antigen receptor-engineered NK cells and of "multifunctional engagers" promoting NK cell activation against tumoral cells further fueled the interest in these lymphocytes (Gardiner, 2019).
NK cell development takes place in the BM and it is largely driven by IL-15 (Kennedy et al, 2000). In the mouse, CD27 single-positive (CD27SP) cells are immature and feature high proliferative and metabolic activity. NK cells further mature through a CD27 and CD11b double-positive (DP) stage to a terminally differentiated CD11b single-positive (CD11bSP) subset, which represents the largest peripheral population and is characterized-unless stimulated-by a quiescent state (Marcais et al, 2014). This is important, as these lymphocytes might otherwise trigger tissue damage or give rise to malignancies (Gianchecchi et al, 2021). These include NK cellderived lymphoma, in which increased levels of the proto-oncogene MYC correlate with poor prognosis (Huang et al, 2014;Wang et al, 2017;Xiong et al, 2020).
MYC is a well-known transcription factor involved in various types of cancers (Dhanasekaran et al, 2022). Even if MYC-binding sites are common throughout the genome, the landscape of its target genes varies based on the investigated cell type (Campbell & White, 2014;Link & Hurlin, 2015). Although the best characterized process controlled by MYC is proliferation, many others are regulated by this transcription factor (Dang et al, 2006;Bretones et al, 2015). In NKT and T cells, its role in development, activation, and proliferation (Bianchi et al, 2006;Dose et al, 2006Dose et al, , 2009Jiang et al, 2010;Nozais et al, 2021) and in orchestrating specific metabolic configurations Marchingo et al, 2020;Saravia et al, 2020) has emerged in the past years.
As for NK cells, exposure to inflammatory cytokines IL-12 and IL-2 or IL-18 has been shown to induce MYC, which in turn increased (G, H) Graphs depict the ratio of Ncr1cre Myc fl/fl over Ncr1cre Myc wt/wt for lymphocytes, total, CD27 SP, DP, and CD11b SP NK cells (CD122 + NK1.1 + CD3/19 − ) in the BM (G) and glycolysis, mitochondrial activity, and effector functions (Loftus et al, 2018;Dong et al, 2019). Stimulation with IL-15, which is crucial for both NK cell homeostasis and effector function, was shown to increase MYC transcription (Cichocki et al, 2009;Ma et al, 2022), which was reported to mediate the up-regulation of NK cell inhibitory receptors (Cichocki et al, 2009) and maintain normal numbers of peripheral NK cells (Dong et al, 2019). Whereas these notions suggest that MYC represents a key factor in NK cells and in mediating the effects of IL-15, we still lack information on its regulation, transcriptional targets, and specific function in these lymphocytes.
We thus set out to investigate MYC expression and the transcriptional landscape regulated by this transcription factor in NK cells. We found that Myc is highly expressed in the more immature BM subsets, whereas its levels are significantly lower in splenic NK cells. In line with this, IL-15 induced and stabilized MYC through distinct mechanisms in both murine and human NK cells. Mice lacking Myc in NK cells presented a substantial reduction of NK cell numbers in the periphery, which was due to a cell-intrinsic defect in developing NK cells. Mechanistically, Myc-deficient NK cells presented a strongly perturbed network of genes coding for ribosome and translation-related proteins, which substantially affected their translation rate and expansion upon IL-15 stimulation, a phenomenon largely recapitulated in human NK cells upon MYC blockade. Together, our results indicate a key role for Myc in regulating translation capacity during IL-15-induced proliferation and the tightly coupled developmental program in NK cells, eventually leading to defective cell numbers and anticancer immunity.

Myc is required for mature NK cell development
We assessed the expression of Myc mRNA by quantitative real-time RT-PCR (qRT-PCR) in BM and splenic NK cells, and splenic NK cells of the DP and CD11b SP subsets ( Fig 1A). Although Myc expression was detected in all tested subpopulations, the highest transcript levels were found in BM NK cells. We therefore tested the expression of Myc at the protein level and found highest Myc expression in CD27SP and DP in the BM, whereas BM-derived CD11b SP and splenic NK cells showed lower levels (Figs 1B and S1A). In line with analyses in human NK cells (Collins et al, 2019), these results suggest that Myc is important for NK cell with high proliferative potential, whereas its reduced levels promote NK cell quiescence in the periphery.
We thus generated Ncr1cre Myc fl/fl mice, in which, Myc expression is deleted in NK cells at the immature stage (Trumpp et al, 2001;Narni-Mancinelli et al, 2011). Deletion was corroborated by qRT-PCR ( Fig S1B). Given that the antibody used for Myc detection can recognize both Myc and n-Myc, generation of Myc-deficient NK cells also allowed us to ensure specificity of the staining for Myc in these lymphocytes ( Fig S1C). Our analyses of Ncr1cre Myc fl/fl mice further showed that the percentage and number of NK cells were significantly reduced as compared with control mice, both in the BM and in the spleen (Figs 1C and D and S1D; a typical gating strategy for NK cells in the BM and spleen is shown in Fig  S1E). Whereas the effect of Myc deficiency on the numbers of CD27 SP NK cells was variable, which is in line with the recent activation of Cre expression in this subset, a prominent decrease in NK cell numbers at the more mature CD11b + maturation stages was observed in the BM and spleen (Figs 1E and F and S1D). Mixed BM chimeras demonstrated that deletion of this transcription factor affected the numbers of CD11b + NK cells in a cell-intrinsic manner (Fig 1G and H). Interestingly, liver conventional NK cells were also markedly decreased in percentage and number, whereas the CD49a + type 1 innate lymphoid cells (ILC1s), in which Ncr1cre-mediated deletion is also expected (Ali et al, 2016;Cuff & Male, 2017;Ducimetiere et al, 2021), did not show a similar decrease ( Fig S1F). Together, these results indicate that Myc is critical for the development of conventional NK cells.

Myc is integral to NK cell proliferation
We thus asked whether the developmental defect of Mycdeficient NK cells was related to a defect in proliferation. We found that the highly dividing CD27 SP NK cells and the DP NK cells from Ncr1cre Myc fl/fl mice presented strongly reduced levels of the proliferation marker Ki67 as compared with the control counterparts in the BM and spleen (Fig 2A and B). Of note, the levels of the antiapoptotic protein Bcl2 were not decreased in these cells ( Fig 2C). Furthermore, we confirmed that reduced percentages of Ki67-positive cells among CD27 SP and DP Mycdeficient NK cells were observed in mixed BM chimeras (Fig 2D and E). This indicates that, in particular, CD27SP NK cells rely on Myc to expand, explaining the decreased numbers of Mycdeficient peripheral NK cells.
Given the importance of IL-15 in driving NK cell development in vivo (Lee et al, 2011;Marcais et al, 2013), we assessed in vitro the response of Myc-deficient NK cells to this cytokine. At an IL-15 concentration not inducing cell division, the survival of these cells was comparable with the one of control NK cells ( Fig 2F). Moreover, in response to IL-15 doses inducing cell division, the percentage of live cells was comparable between the two genotypes, whereas Myc-deficient NK cells exhibited severely compromised expansion, hinting to an important proliferative defect ( Fig 2F). Of note, a similar result was obtained also in response to IL-2, a cytokine sharing two receptor subunits with IL-15 ( Fig 2F).

Myc is regulated by multiple IL-15-triggered signaling pathways
Having observed that Myc-deficient NK cells exhibit an expansion defect in response to IL-15 and IL-2 stimulations, we measured the expression of the two shared receptor chains CD122 (IL2Rβ/IL15Rβ) and CD132 (common γ chain). CD122 showed a moderately reduced expression on Myc-deficient NK cells, whereas we did not detect differences in CD132 levels (Fig S2A and B). Next, we thought to assess the phosphorylation of signal transducer and activator of transcription 5 (STAT5), extracellular signal-regulated kinase (ERK), and S6, key signaling pathways induced by IL-15 in these cells. Multiple cascades were affected in the more mature Myc-deficient NK cell subsets ( Fig 3A). Being the earliest deleted subpopulation and the one with the strongest proliferation defect (Fig 2A, B, D, and E), we were however particularly interested in the CD27 SP subset, which only presented a reduction in ERK signaling noticeable already ex vivo ( Fig 3A). Immunoblot analysis of splenic Myc-deficient NK cells did not reveal major differences in the expression of total STAT5, S6 or ERK2 ( Fig S2C).
As previous results showed that an even stronger ERK signaling defect in Shp-2-deficient NK cells was not associated to an impaired development (Niogret et al, 2019), we decided to investigate if instead IL-15 was regulating Myc. In line with previous findings, we observed that IL-15 stimulation markedly increased Myc transcript levels ( Fig 3B) (Cichocki et al, 2009;Mishra et al, 2012;Gotthardt et al, 2016;Pinz et al, 2016;Villarino et al, 2017). The inhibitors tofacitinib (pan JAK inhibitor), PD98059 (mitogen-activated protein kinase kinase inhibitor), and Torin2 were used to block the JAK, ERK, and mammalian target of rapamycin (mTOR) pathways, respectively ( Fig  S2D) (Pende et al, 2004;Meyer & Levine, 2014). Whereas tofacitinib and Torin2 affected Myc mRNA induction, concomitant treatment almost prevented it ( Fig 3B). Stimulation of splenic NK cells also increased Myc protein levels ( Fig 3C). To distinguish the contribution of the posttranslational regulation of Myc from its transcriptional  (A) Volcano plot displaying differential expressed genes; the vertical axis (y-axis) displays the significance as the log 10 (P-value), and the horizontal axis (x-axis) displays effect size as the log 2 fold change value (FDR = 0.01; logFC threshold = 0.5; genes shown are significant for the three following statistical methods; DESeq2 [Wald test], edgeR, limma-trend). (B) Hierarchical clustering of the 138 significantly differentially expressed genes. (C) The table shows the top significant KEGG gene sets (q ≤ 0.0006 with fGSEA, fisher, and GSVA methods-indicated under Stars), their fold-change, and the number of genes in each. Enrichment plot for ribosome KEGG gene set is shown. induction, we treated IL-15-stimulated NK cells with cycloheximide to prevent further protein synthesis and studied the effects of mTOR and ERK, two pathways implicated in Myc stabilization (Sears et al, 2000;Gustafson & Weiss, 2010), on Myc levels ( Fig 3D). We observed that Myc degradation was accelerated by the blockade of mTOR, whereas ERK pathway inhibition exerted negligible effects. These results indicate that, in IL-15-treated NK cells, JAK-STAT and mTOR are mostly responsible for Myc induction, whereas mTOR for Myc stability.

Myc regulates genes related to translation
To delineate the role of Myc in NK cell expansion and development, we performed RNA sequencing (RNA-seq) analyses. Being the earliest deleted subset and presenting the strongest proliferation defect (Fig 2A, B, D, and E), we selected Myc-deficient and -sufficient CD27SP NK cells harvested from the BM. Paired differential expression analyses identified 69 genes significantly down-regulated and 69 genes significantly up-regulated in the absence of Myc ( Fig  4A). We next performed a hierarchical clustering of these genes ( Fig  4B). Notably, 36 of the down-regulated genes coded for ribosomal subunits, translation elongation factors, and proteins involved in ribosomal formation ( Fig 4B; indicated in bold). Among the pathways revealed by Gene set enrichment analysis with highest significance (Kanehisa & Goto, 2000), the most deregulated referred to "Ribosome" (Fig 4C). Of note, this finding is in line with work by Tang et al (2021), who demonstrated a defect in translation-related genes in Myc-deficient bulk splenic NK cells. In keeping with translationrelated genes, we also identified Slc7a5, a Myc target involved in amino acid transport, among the down-regulated genes in Mycdeficient immature NK cells ( Fig 4B) (Marchingo et al, 2020). In line with our flow cytometric analyses, we observed down-regulated expression of the gene coding for Ki67 (Mki67), which has also been linked to nucleologenesis ( Fig 4B) (Schmidt et al, 2003;Rahmanzadeh et al, 2007). Instead, among these significantly altered genes, we did not observe canonical Myc-responsive genes involved in cell cycle. We therefore specifically looked for their expression, revealing that part of them was tendentially-but not significantly-altered, whereas others were not ( Fig S3A) (Robson et al, 2011;Mishra et al, 2012;Bretones et al, 2015). Similar results were observed for other pathways classically associated to Myc, including oxidative phosphorylation, glycolysis, glutamine, and polyamine ( Fig S3B-E) (Bult et al, 2019). These data reveal a prominent role for Myc in ribosomagenesis and translation more broadly in NK cells, possibly overriding its most established function in directly regulating cell cycle progression.

Sustained translational capacity is essential for NK cell expansion
We therefore tested the impact of Myc absence on the translation rate of NK cells. To study this process, we assessed the incorporation in newly synthesized proteins of 35S-radiolabeled methionine and cysteine by splenic NK cells exposed for an overnight period to limited or high IL-15 concentration. When acutely stimulated, Myc-deficient NK cells exhibited defective translation ( Fig  5A). We therefore verified selected RNA-seq results by Western blot and confirmed that multiple ribosomal proteins were expressed to a lower extent in Myc-deficient NK cells (Fig 5B). This led us to investigate the effects of low doses cycloheximide, aimed at partially inhibiting translation, on the expansion of WT NK cells. Recapitulating the phenotype of Myc-deficient NK cells, low doses of CHX did not affect NK cell survival at low IL-15, but impaired NK cell expansion in response to high IL-15 concentration (Fig 5C and D). Furthermore, cycloheximide and Myc deficiency limited the raise in NK cell size and granularity observed upon exposure to acute IL-15 stimulation (Fig 5C and E). Taken together, our data underline how a partial defect in translation phenocopies the defects caused by Myc deficiency.

IL-15 regulates MYC and ribosomal proteins also in human NK cells
In murine NK cells, we observed that IL-15 up-regulated Myc through several mechanisms. We thus aimed to investigate whether this phenomenon was true in human NK cells. For this, we relied on DERL7, a human lymphoma cell line featuring NK cell lineage characteristics (Di Noto et al, 2001;Dufva et al, 2018). Indeed, IL-15 stimulation markedly increased MYC transcript and protein levels ( Fig 6A and B). Inhibition of the JAK and mTOR pathways and, in particular, concomitant treatment with tofacitinib and Torin2 affected MYC mRNA induction ( Fig 6A). Similar to mouse NK cells, in DERL7 cells, MYC degradation was accelerated by the blockade of mTOR, whereas ERK pathway inhibition did not affect MYC stability ( Fig 6C). These data highlight strong analogies between the regulation of MYC in human and murine NK cells.
To evaluate whether MYC promotes ribosome formation and translation also in human cells, we isolated NK cells from healthy donors' blood and assessed the effects of the specific MYC inhibitor 10058-F4 upon stimulation with 0.25 and 25 ng/ml human IL-15 in vitro. Human NK cells did proliferate significantly already at the low dose of IL-15 and MYC inhibition limited their expansion most prominently in response to high IL-15 stimulation (Figs 6D and E and S4A). Notably, similar effects were achieved by treating human NK cells with low-dose CHX, aimed at partially reducing translation (Figs 6D and E and S4A). Whereas the difference in size and granularity induced on human NK cells going from 0.25 to 25 ng/ml IL-15 was marginal, 10058-F4 and CHX also showed a tendency to reduce NK cell size and granularity (Figs 6F and G and S4A). Notably, 10058-F4 treatment decreased the expression of the tested ribosomal proteins (Fig 6H) induced by-in particular-25 ng/ml of IL-15. These results underline strong parallels between mouse and human NK cells with respect to MYC regulation and its effects on the translation machinery.
We next asked whether ribosomal proteins were affected at the transcriptional level also in primary human NK cells after MYC inhibition. To this aim, we treated primary NK cells with IL-15 in presence or absence of 10058-F4 and checked the mRNA levels of selected ribosomal genes. Transcripts of these genes were however not reduced (Fig S4B). We thus wondered if this could be attributed to the chemical inhibitor used in human NK cells. To assess this, we compared ribosomal gene transcripts in murine control NK cells treated with a Myc inhibitor and in Myc-deficient NK cells. 10058-F4 was not only less effective than Myc ablation to reduce the levels of these transcripts, but also showed transient effects (Fig S4C), indicating that the inhibitor works inefficiently or on some aspects of MYC function. In spite of this, MYC inhibition culminated in affecting the translation machinery at the protein level (Fig 6H), suggesting that MYC regulates translation in a robust and multilevel fashion.

Deletion of Myc impacts on NK cell immunity
We finally asked whether the perturbations caused by the absence of Myc on NK cell homeostasis influenced immunity. Given the importance of inhibitory and activating receptors for NK cell education and function, we first assessed the receptor repertoire on NK cells from Ncr1cre Myc fl/fl and Ncr1cre Myc wt/wt mice. Although a decrease in the levels of KLRG1 on the effector subsets was noticeable, most other differences were moderate or absent (Fig 7A). A second parameter essential for NK cells' function is their capacity to produce effector cytokines and cytotoxic mediators. As IL-15 is not only a crucial cytokine for NK cell ontogeny and homeostasis, but also for NK cell effector capacity, we measured the production of cytotoxic or inflammatory mediators by Myc-deficient NK cells before or after exposure to IL-15. Although the production of granzyme (Gzm) A was high irrespective of the genotype or the stimulatory condition in DP and CD11b SP NK cells, we observed a substantial defect in the production of GzmB and IFN-γ in DP and CD11b SP Myc-deficient NK cells upon IL-15 stimulation (Fig 7B). These results are in good agreement with the curtailed response to IL-15 that we observed in these NK cell subsets (Fig 3A).
To understand if these differences resulted in a decreased killing capacity by Myc-deficient NK cells, we examined in vitro their cytotoxic activity. Myc-deficient NK cells were able to significantly eliminate RMA-S (lacking MHCI) target cells, despite to a lower extent when compared with control NK cells, a trend that was observed also in the aspecific toxicity affecting control RMA cells (expressing MHCI) (Fig 7C).
Finally, we studied how deletion of Myc affects immunity in vivo. We thus injected Ncr1cre Myc fl/fl and Ncr1cre Myc wt/wt mice i.v. with B16 cells, a tumor model for metastatic melanoma. Notably, based on the number of nodules found 10 d later in the lung, Ncr1cre Myc fl/fl mice eliminated significantly less melanoma cells (Fig 7D). Taken together, these data highlight how Ncr1cre Myc fl/fl mice present a substantially weakened antitumor immunity, reflecting a moderate NK cell-intrinsic killing defect together with a substantial loss of mature peripheral NK cells.

Discussion
Here, we report that Myc is critical for NK cell development and, thus, for NK cell-mediated immunity. With regard to NK cell effector functions, we observed that Myc was needed for the induction of IFN-γ and GzmB by IL-15, in line with findings downstream of IL-2/IL-12 stimulation (Loftus et al, 2018). Indeed, we found that Myc was required for the IL-15-dependent activation of ERK, S6, and STAT5 in effector NK cells. Despite this, neither the levels of granzyme A nor the killing on a per cell basis were dramatically affected by Myc deletion. Furthermore, with the exception of KLRG1, peripheral Mycdeficient NK cells exhibited mild alterations in the expression of inhibitory receptors. These findings did not support a prime role for Myc in controlling cytotoxicity and inhibitory receptor induction in our system. In Ncr1cre Myc fl/fl mice, peripheral mature NK cell numbers were decreased by roughly 70%, largely explaining the defective NK cell-mediated anticancer immunity.
We thus demonstrated a major role for Myc in developing NK cells, in which it is highly expressed. In line with these findings and with previous publications on the role of Myc in IL-15-dependent T cell homeostasis and development (Bianchi et al, 2006;Jiang et al, 2010), we showed that Myc is required for NK cell expansion upon IL-15 stimulation. As IL-15 normally induced multiple signaling pathways in CD27SP Myc-deficient NK cells, we asked whether Myc was implicated downstream of STAT5, mTOR, and ERK. Besides corroborating the role of STAT5 (Gotthardt et al, 2016;Villarino et al, 2017), we identified a contribution by mTOR in regulating MYC transcript levels. Furthermore, we found that mTOR is crucial to limit MYC degradation downstream of IL-15. Altogether, these results not only demonstrate that MYC is a central node downstream of IL-15, but also identify the multilevel crosstalk between MYC and IL-15 in murine and human NK cells.
In IL-2/IL-12-stimulated NK cells, Myc has been shown to promote glycolysis and mitochondrial function (Loftus et al, 2018). While we confirm a trend towards decreased expression of genes involved in these pathways in Myc-deficient NK cells, we found that most of the genes significantly regulated by Myc are related to translation, a fundamental cellular process, which might-in turn-broadly affect cellular metabolism. The capacity of Myc to transactivate ribosomal genes has so far been mostly studied in cancer and pluripotent cells (Poortinga et al, 2004;van Riggelen et al, 2010;Seitz et al, 2011;Lee et al, 2012;Scognamiglio et al, 2016;Thomas et al, 2019;Destefanis et al, 2020;Popay et al, 2021). Supporting a role for Myc in basic protein metabolism, a recent study on activated T cells highlighted how Myc controls the expression of amino acid transporters including Slc7a5, whose deletion largely phenocopied the effects of Myc ablation (Marchingo et al, 2020). In NK cells, although we found that Myc maintains Slc7a5 transcription, Slc7a5 is also required to sustain Myc levels (Loftus et al, 2018), revealing a feed-forward loop concurring to regulate translation. Notably, emerging data show that translation-related genes represent a key but often overlooked signature controlled by Myc in multiple primary immune cells and future work will help define the important roles that such a basic mechanism might play in these cells (Perez-Olivares et al, 2018;Marchingo et al, 2020;Saravia et al, 2020;Tang et al, 2021).
Given the impossibility to genetically rescue the multiple defects of the translational module observed in Myc-deficient NK cells, we explored the effects of partially inhibiting translation in NK cells. Notably, this approach recapitulated the expansion, size, and granularity phenotype observed in both murine and human NK  However, we also observed that using 10058-F4 to block MYC activity in vitro did not fully mirror the phenotype of Myc-knockout NK cells at the transcriptional level. This suggests that this inhibitor acts either less effectively or differently from genetic deletion. Yet, we conclude that MYC regulates the translational machinery through multiple pathways that remain to be elucidated in greater detail. Importantly, we also identified cycloheximide doses which affected proliferating NK cells but not resting NK cells' viability, revealing an interesting therapeutic window.
Our results suggest that drugs blocking translation, IL-15, JAK or mTOR downstream of this cytokine might concur in the treatment of lymphoid tumors of NK cell origin depending on MYC (Devlin et al, 2016;Schwartz et al, 2017;Laham-Karam et al, 2020). Our data also indicate that the low levels of Myc found in unstimulated peripheral NK cells limit their protein synthesis rate, thereby warranting a status of quiescence, which is important to avoid overzealous responses or a state of terminal differentiation (Castro et al, 2018).
Although the contribution of Ncr1 + ILCs, ILC1s in particular, to the control of B16 melanoma cannot be formally excluded, our observation that liver ILC1s are not reduced and the notion that NK cells are key to control this metastatic cancer model supports the relevance of Myc in conventional NK cells for maintaining the antitumor response (Victorino et al, 2015;Lopez-Soto et al, 2017;Howard et al, 2021). This observation is of great importance as in the work by Tang et al, no significant difference was observed when subcutaneously engrafting B16 melanoma cells in Ncr1cre Myc fl/fl mice as compared with the controls (Tang et al, 2021).
Whereas our results suggest that combining NK cell-based immunotherapeutic approaches with strategies aimed at targeting MYC might be inappropriate, data by others suggest that targeting MYC in lymphoma cells exerts beneficial effects on NK cell responses (Swaminathan et al, 2020). Given these complementary observations, it will be important to weight the effects of Myc inhibitors in preclinical cancer models and evaluate different regimens (Han et al, 2019). The multilevel interplay between MYC and IL-15 and the role of this transcription factor in anticancer NK and-most likely-chimeric antigen receptor-NK cells' immunity reinforce the rationale for developing and employing targeted immunostimulatory approaches, such as multifunctional engagers carrying IL-15 (Lee et al, 2011;Vallera et al, 2016;Liu et al, 2018;Gardiner, 2019;Zakiryanova et al, 2019).

Human primary NK cell isolation and culture
Human peripheral blood mononuclear cells (PBMCs) from healthy donors were prepared form purchased buffy coats (Servizio Trasfusionale della Svizzera italiana; Lugano). Buffy coats were diluted 1:2 by adding an equal amount of prewarmed RPMI 1640 (Gibco) and layered over Ficoll (Cytiva). Cells were span down at 800g for 30 min, at 20°C, without break. Interphase with PBMCs was harvested and diluted with a prewarmed medium followed by spinning down at 450g for 15 min at 20°C and washing with the above medium containing 0.5% human serum. NK cells were enriched from PBMCs by EasySep Human NK Cell Isolation Kit (Stem Cell Technologies) and cultured at 37°C in 5% CO 2 in RPMI 1640 (Gibco) supplemented with 10% FBS (Gibco), 50 U/ml penicillin, 50 μg/ml streptomycin (Gibco), 1X MEM NEAA (Gibco), 1 mM sodium pyruvate (Gibco) in the presence of the indicated amounts of human recombinant IL-15 (Peprotech).

Mice
Ncr1cre Myc fl/fl and Ncr1cre Myc wt/wt mice were generated by crossing Ncr1cre and Myc fl/fl mice which were previously described (Trumpp et al, 2001;Narni-Mancinelli et al 2011). Ncr1cre Myc wt/wt mice were also crossed onto a CD45.1 congenic background. The mice were backcrossed multiple times on a C57BL/6 background. These mice and C57BL/6 mice were crossed and/or housed at the animal facility of the University of Lausanne and of the Institute for Research in Biomedicine of Bellinzona. Sex-and age-matched 6-14-wk-old mice were used in the experiments. All animal experimental protocols were approved by the veterinary office regulations of the states of Vaud and of Ticino, Switzerland, and all methods were performed in accordance with the Swiss guidelines and regulations.

Isolation of liver immune cells
Livers were perfused with PBS 10 IU/ml heparin (Sigma-Aldrich) and mechanically dissociated by gently pressing through 40 μM nylon mesh filters. The cell suspension was centrifuged at 50g for 5 min to remove the parenchymal cells. The supernatants were centrifuged for 5 min at 500g and the pellet was washed twice with RPMI 1640 (Gibco) medium. The pellet was then resuspended in 40% Percoll (Cytiva) in HBSS (Gibco) and gently overlaid onto 70% Percoll and centrifuged at 800g for 25 min without a break. The cells were then collected at the interface and washed with PBS and used for flowcytometry analysis.   RPMI-1640 (Gibco) supplemented with 10% FBS (Gibco), 2 mM L-glutamine (Gibco), 50 U/ml penicillin, 50 μg/ml streptomycin (Gibco) and cultured at 37°C in 5% CO2. Then, 10 5 cells were intravenously injected in Ncr1cre Myc wt/wt mice and Ncr1cre Myc fl/fl mice. The mice were euthanized at day 10 after injection and the number of lung surface nodules was counted.

Mixed BM chimeras
Mixed BM chimeras were generated as previously described (Niogret et al, 2019), except that donor BM cells were obtained from sex-matched Ncr1cre Myc fl/fl mice and Ncr1cre Myc wt/wt mice with the CD45.1 congenic marker by flushing tibias and femurs and mixing them in a 1:1 ratio before injection. The chimeras were used for the respective experiments 6-8 wk post-BM reconstitution.
Streptavidin-conjugated fluorophores were purchased from BioLegend or eBioscience. Stainings were performed with appropriate combinations of fluorophores. Propidium iodide (0.05 μg/ml from Immunochemistry) was used as a live-dead marker. Events were collected on a FACS Canto, Fortessa (Becton Dickinson) and data were analyzed using FlowJo software (Tree Star Inc.).

Intracellular flow cytometry stainings
For intracellular staining of Ki67 (SolA15, 1:100; eBioscience), murine NK cells were first surface stained, fixed, and permeabilized using the FoxP3 transcription factor staining buffer set of eBioscience (00-5523-00) according to the manufacturer's recommendations. For assessment of Myc levels, intracellular staining was performed using Foxp3/Transcription Factor Staining Buffer Set (eBioscience) according to the manufacturer's protocol. Cells were stained for 1 h in a permeabilization buffer with an Myc antibody (Myc/n-Myc D3N8F rabbit mAb PE antibody, #35876, 1:50; Cell Signaling). Myc expression (expressed as geometric MFI) was measured by flow cytometry and the background of Myc-knockout NK cells was subtracted. For Bcl2 staining, BM cells were fixed, permeabilized as above, and stained for 1 h in the permeabilization buffer with Bcl2 antibody (Bcl2 10C4, 1:100; Invitrogen).
To assess the Myc stability, enriched mouse splenic NK cells or human DERL7 NK cell line were stimulated 2 or 3 h, respectively, with recombinant mouse or human IL-15. The cells were then treated with CHX (20 μg/ml) alone or in combination with Torin2 (250 nM from Selleckchem), PD98059 (10 μM) or both for 70 min. The Myc content was then assessed by intracellular staining as mentioned above.

Analyses of effector molecule production
For intracellular IFN-γ, GzmA, and GzmB detection, splenocytes from control and Ncr1cre Myc fl/fl mice were left unstimulated or stimulated for a total of 4.5 h at 37°C with 50 ng/ml IL-15 and Brefeldin A (10 μg/ml Enzo Life Science) was added after the first 2 h. The cells were fixed with fixation/permeabilization buffer (eBioscience) for intracellular staining for IFN-γ (F3 IGH48; eBioscience) and granzyme A (3G8.5; eBioscience). Granzyme B (NGZB; eBioscience) staining was performed after fixation (25 min at 22-28°C) and permeabilization with the Foxp3 Transcription Factor Staining Buffer Set (#00-5523-00; eBioscience). Intracellular staining was then performed for 30 min at 22-28°C (with the antibody diluted in the permeabilization buffer).

Mouse NK cell isolation and culture
For experiments requiring enriched NK cells, splenic NK cells were isolated using the EasySep mouse NK cell isolation kit (Stemcell Technologies) according to the manufacturer's recommendations. NK cell enrichment was confirmed by flow cytometry. As indicated, NK cells were further FACS-sorted for selected experimental settings and when derived from BM. NK cells were then cultured in RPMI 1640 (Gibco) supplemented with 10% FBS (Gibco), 50 U/ml penicillin, 50 μg/ml streptomycin (Gibco), 1 mM sodium pyruvate

Metabolic labeling
Mouse NK cells (purity > 94%) were washed with prewarmed (PBS; Gibco) and then pulse labeled with 0.05 mCi of [35S]-methionine/ cysteine in DMEM (Gibco), 50 mM Hepes (Gibco), and 2 mM L-glutamine (Gibco) for 20 min at 37°C. The labeling medium was removed and the cells were washed with PBS. The cells were lysed with RIPA buffer (1% Triton X-100, 0.1% SDS, 0.5% sodium deoxycholate in HBS, pH 7.4) containing 20 mM N-Ethylmaleimide and protease inhibitor cocktail (1 mM PMSF, chymostatin, leupeptin, antipain, and pepstatin, 10 μg/ml each) for 20 min on ice. Lysates were subjected to centrifugation at 4°C and 10,000g for 10 min and postnuclear supernatants were collected. Samples were subjected to SDS-PAGE, and proteins were transferred to PVDF membranes using the Trans-Blot Turbo Transfer System (Bio-Rad). PVDF membranes were exposed to autoradiography films (GEHealthcare) and were scanned with the Typhoon FLA 9500 (Software Version 1.0).

Cell labeling and in vitro expansion analyses
To study murine NK cell expansion, enriched NK cells were labeled with 5 μM cell trace violet (CTV; Life Technologies) in PBS 2% FBS at 37°C for 20 min. Cells were seeded at a density of 0.6-2 × 10 6 /ml in the presence of the indicated concentration of recombinant mouse IL-15. In the experiments where cycloheximide (CHX; Sigma-Aldrich) was employed, after an overnight period, 0.01 or 0.05 μg/ml of CHX was added and the cells were further cultured for 3 d.
To study human NK cell expansion, enriched NK cells were labeled with 5 μM CTV, seeded at a density of 0.5 × 10 6 /ml and stimulated with 0.5 or 25 ng/ml recombinant human IL-15 (PeproTech). After an overnight period, 0.04 μg/ml of CHX, 20 μM MYC inhibitor 10058-F4 (Sigma-Aldrich) was added, and the cells further cultured for 3 d. The cells were stained and analyzed by flow cytometry.
Quantitative RT-PCR analyses BM and splenic NK cells were FACS-sorted as indicated in the text. Sorts were performed using FACSAria III, BD Biosciences. In other experiments, isolated WT NK cells were pretreated (15-30 min) or not with Torin2 (250 nM from Selleckchem), PD98059 (10 μM from Adipogen), tofacitinib (90 nM from GSK) or combinations thereof and stimulated with 50 ng/ml of recombinant mouse IL-15 for 2 h. Then, total RNA was extracted using the TRIzol reagent (Life Technologies) according to the manufacturer's instructions. Annealing with random primers (Life technologies) was performed at 70°C for 5 min, followed by retrotranscription to cDNA with M-MLV RT, RNase H(−) point mutant (Promega) in the presence of buffer, nucleotides, and RNasin Plus RNase Inhibitor (Promega). The reaction was then incubated at 40°C for 10 min, 45°C for 50 min, and 70°C for 15 min. cDNAs obtained were diluted and quantitative PCR was performed using the PerfeCTa SYBR Green FastMix (Quantabio) on a QuantStudio 3 Real-Time PCR System machine (Applied Biosystems). Standard cycling was used (45 cycles of 95, 60, and 72°C steps of 10 s each). Negative control reactions were cycled alongside test samples to ensure the absence of contaminating genomic DNA. The expression of Myc gene was determined relative to the abundance of the housekeeping gene RNA polymerase II Subunit A (Polr2a) and/or 18 s. Data were analyzed and transcript abundance and SD calculated using the Thermo Fisher Connect software.

Library preparation and sequencing
Total RNA was extracted from sorted cells with an RNeasy plus mini kit (QIAGEN) according to the manufacturer's instructions.
The quantity and quality of the isolated RNA was determined with a Qubit (1.0) Fluorometer (Life Technologies) and a 4200 TapeStation System (Agilent). The libraries were prepared following the Takara SMART-Seq v4 Ultra Low Input RNA Kit, with PolyA selection. This protocol relies on the template switching activity of reverse transcriptases to enrich for full-length cDNAs and to add defined PCR adapters directly to both ends of the first-strand cDNA. Briefly, total RNA samples (0.25-10 ng) were reverse transcribed using random priming into double-stranded cDNA in the presence of a template switch oligo. When the SMARTScribe Reverse Transcriptase reaches the 59 end of the RNA, its terminal transferase activity adds a few nucleotides to the 39 end of the cDNA. The SMART-Seq Oligonucleotide base-pairs with the non-templated nucleotide stretch, creating an extended template to allow the reverse transcriptase to continue replicating to the end of the oligonucleotide. The SMART-Seq primer and oligo serve as universal priming sites for cDNA amplification by PCR. 0.5 ng of cDNA from each sample were tagmented and amplified using Illumina Nextera XT kit. The resulting libraries were pooled, double-sided size selected (0.5× followed by 0.8× ratio using Beckman Ampure XP beads), and quantified using an Agilent 4200 TapeStation System.
For the CD27SP NK cell isolated from the BM, after quantification, libraries were prepared for loading accordingly to the NovaSeq workflow with the NovaSeq6000 Reagent Kit (Illumina). Cluster generation and sequencing were performed on a NovaSeq6000 System with a run configuration of single end 100 bp with a depth of around 20 Mio reads per sample. FASTQ files were processed using Nextflow and the nfcore/RNAseq pipeline with default values (Di Tommaso et al, 2017;Ewels et al, 2020). Quantification of transcript expression on the gene level was computed using Salmon and the GRCh38 reference transcriptome (Patro et al, 2017).

RNA-seq analyses
For gene-level testing, statistical significance was assessed using three independent statistical methods to be more robust: DESeq2 (Wald test), edgeR (QLF test), and limma-trend (Robinson et al, 2010;Love et al, 2014;Ritchie et al, 2015). The maximum q-value of the three methods was taken as aggregate q-value, which corresponds to taking the intersection of significant genes from all three tests. Statistical testing and visualization (heatmap, volcano plot) have been performed using the Omics Playground version v2.7.8 (Akhmedov et al, 2020). Gene set enrichment analysis was performed using the KEGG gene sets (Kanehisa & Goto, 2000) and a combination of three statistical tests: fGSEA (Korotkevich, 2021 Preprint), (Fisher, 1922), and GSVA (Hanzelmann et al, 2013).

Data Availability
RNA-seq datasets are available in the GEO repository under the number GSE220930. The data that support the findings of this study are available from the corresponding author upon reasonable request.
A Trumpp: resources. G Guarda: conceptualization, data curation, supervision, funding acquisition, project administration, and writing-original draft, review, and editing.